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Abstract 

We extend the multifractal analysis of the statistics of critical wave functions 
in quantum Hall systems by calculating numerically the correlations of local 
amplitudes corresponding to eigenstates at two different energies. Our results 
confirm multifractal scaling relations which are different from those occurring 
in conventional critical phenomena. The critical exponent corresponding to 
the typical amplitude, ao ~ 2.28, gives an almost complete characterization of 
the critical behavior of eigenstates, including correlations. Our results support 
the interpretation of the local density of states being an order parameter of 
the Anderson transition. 
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I. INTRODUCTION 



Two-dimensional independent electrons in the presence of static disorder and a strong 
magnetic field undergo disorder induced localization-delocalization (LD) transitions when 
the Fermi energy crosses critical energies - the Landau energies. These LD transitions are 
believed to be responsible for the integer quantum Hall effect In finite size systems 

the localization length £ of the electronic states is larger than the system size L for a certain 
energy range, AE, around the critical energies. These states are called critical states. In the 
thermodynamic limit AE oc L~ x l v where v is the critical exponent of £. An obvious task for 
a theory of the LD transition is to yield the statistics and scaling behavior of eigenstates in 
the critical regime. After pioneering works by Wegner fU and Aoki @] it became clear that 
the critical wave functions have a multifractal structure (for a review see |J and references 
therein). The entire distribution of local amplitudes and its scaling behavior is encoded in 
the multifractal f{ct) spectrum. The distribution is broad on all length scales and close to a 
log-normal distribution. The most important quantity is the maximum position, ato, of f(a). 
It describes the scaling behavior of the geometric mean of what serves as a typical amplitude 
of a critical wave function. The f(a) spectrum has been interpreted as a spectrum of critical 
exponents related to the order parameter field of the LD transition ||. To be consistent 
with such an interpretation f(a) has to share several features with scaling exponents in 
conventional critical phenomena. Firstly, it has to be universal, i.e. independent of the 
disorder configuration and the microscopic details of the electrons state. In previous studies 
this was confirmed ||. Secondly, correlations of the order parameter field have to be related 
to f(a) by appropriate scaling relations. A systematic investigation of this topic is presented 
here (Related conformal scaling relations were recently investigated by Dohmen et al. @). 
Our numerical data are consistent with universal scaling relations between f{ot) and scaling 
exponents of the energetic and spatial correlations of critical wave functions. 

The article is organized as follows. In Sec. H we explicitly demonstrate that the mul- 
tifractal spectrum describes the distribution function of local amplitudes. In Sec. [I I| the 
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spatial correlations of local amplitudes (and powers thereof) of wave functions at a fixed 
energy in the critical regime are investigated. We study the scaling with respect to the 
size of boxes which the local amplitude is averaged over, and with respect to the distance 
between correlated amplitudes. We find that scaling relations are fulfilled which can also be 
obtained by heuristic arguments. In Sec. [TV] we study the correlator of two local amplitudes 
corresponding to critical eigenstates at energies separated by u. We find that a characteristic 
length scale (introduced by Chalker f|||) serves as a cutoff length for multifractal corre- 
lations. Remarkably, two local amplitudes corresponding to eigenstates separated in energy 
by uj are correlated in the same way as two amplitudes corresponding to one eigenstate, 
provided their spatial distance is less than L w . By transforming from the wave function 
amplitudes to the local density of states our results show that the local density of states has 
several features in common with order parameter fields. The conclusions are summarized in 
Sec. 

II. DISTRIBUTION OF LOCAL AMPLITUDES 

In this section we first review the multifractal analysis of the distribution of local am- 
plitudes of critical eigenstates. Secondly, we confirm that the distribution is encoded in the 
multifractal f(a) spectrum by a direct comparison of the numerically obtained histogram of 
local amplitudes with the theoretical distribution function that follows from the numerically 
obtained f(a) spectrum. 

In 1983 Aoki [|J gave a nice argument for the multifractal behavior of critical wave 
functions (although at that time the phrase 'multifractality' was not yet common). His 
argument goes as follows. Consider the inverse participation number defined by 

V = J d d r\^(r)\ 4 (1) 
n 

where Q denotes a (i-dimensional region with linear size L. If the wave function ip(r) is 
uniformly distributed - like in a metallic phase - then V oc L~ d and the participation 
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ratio p = (VL d )^ 1 is constant. In the localized regime V ~ and p vanishes in the 
thermodynamic limit. At the transition point where the wave function is extended the 
participation ratio still has to vanish in the thermodynamic limit if the LD phenomenon is 
similar to a second order phase transition. Consequently, V scales with a power d* < d. 
Wegner had already calculated the whole spectrum of exponents for generalized inverse 
participation numbers within the non-linear sigma-model |§. This spectrum was interpreted 
as a multifractal spectrum by Castellani and Peliti [|K|] . After extensive numerical work the 



following description of the statistics of critical wave functions ip(r) is now established: 
Consider the box-probability 

P(U) := / |^(r)| 2 (2) 

J box 

of some box with linear size L b , normalized to the total volume L d , P(L) = 1. At the LD 
transition the corresponding distribution function ir(P] L b /L) gives rise to power law scaling 
for the moments, 

([P^K^, (3) 

where d + r(g) is a non-linear function of q. This non-linearity is a direct consequence of 
Aoki's observation that d + r(2) = d + d* ^ d + d. The brackets (■■■) L in Eq. denote 
the average over disorder configurations. In practice it turns out that, to a good accuracy 
this average can often be substituted by the spatial average over one wave function for a 
given configuration. Within numerical accuracy the resulting spectra are identical. This 



"universality" is expected to be precise in the thermodynamic limit (see also [11]). 

The corresponding (universal) distribution function can be described in terms of a single- 
humped, positive function f(a) called the multifractal spectrum of the wave function, 

tt(P; L b /L) dP oc (L b /L) d ~ f{a) da , (4) 

where a := InP/ ]n(Lf,/L); /(a) is related to r(q) by a Legendre transformation 

f(a(q)) = a(q)q- r(q) , a(q) = dr(q)/dq . (5) 



Thus, the statistics of critical eigenstates is encoded in f(a) or equivalently in r(q). Because 
of the horse-shoe shape of f(a) it can be approximated by a parabola, 

f(a) ^ d- (a - a ) 2 /(4(a - d)) . (6) 

This parabolic approximation (PA) contains «o as the only parameter besides d. This is 
due to the assumed validity of the PA at least up to \q\ < 1. Equation (g) corresponds to a 
log-normal distribution centered around the typical value P typ '■— exp < InP >oc (Lb/L) a ° 
with log- variance proportional to — d. A simple one-parameter approximation for f(a) 
which takes into account that the support [a(oo), a(— oo)] of f(a) is finite, is the semi-elliptic 
approximation (SEA) 



f(a) « d, 



(a - «o) 2 ^ 



\ «o — d 2 

To demonstrate that the distribution of local amplitudes of critical eigenstates is encoded 
in f(a) we present numerical results for a quantum Hall system (QHS). 

The wave functions are calculated for the model of independent (spinless) electrons sub- 
ject to strong magnetic field and disorder. The disorder was implemented by a set of 5- 
impurities with random positions and random strengths symmetric around zero. To avoid 
that some of the (degenerate) wave functions of the pure system have zeroes precisely at the 
position of the point scatterers (and thus are not affected by disorder) the number of point 
scatterers was taken to be larger then the number of flux-quanta penetrating the system's 
area L 2 (cf. ||12|| ). The microscopic length scale of the problem is the magnetic length Is 
defined by the size of a cell penetrated by a single flux-quantum, 2nl 2 B . We worked out the 
representing Hamiltonian matrix in the Landau representation which is convenient for peri- 
odic boundary conditions in one direction (say y— direction). To account for a finite system 
size in x— direction we adopted the Landau counting procedure which results in a matrix 
of dimension N = L 2 /(2nl 2 3 ), after projecting to the lowest Landau band. The matrix has 
band structure with a bandwidth of order y/N and allows for diagonalizing systems of lin- 
ear size L about 2001b with the aid of usual workstations. The diagonalization yields the 



eigenvalues and eigenstates for any desired energy window within the lowest Landau band. 
The determination of the range AE of critical states was based on a previous analysis of 



Thouless numbers within the same model 13 



In Fig. [I] the squared amplitudes of a wave function from the center of the Landau 
band are shown together with the f(u) spectrum calculated from these amplitudes. The 
corresponding histogram of the logarithm of amplitudes (measured on a box of size 4Z^) is 
displayed in Fig. || together with the distribution function calculated from the f(a) spectrum 
using Eq. (J|). These figures demonstrate that the distribution of amplitudes is (i) encoded 
in the f(a) spectrum and (ii) is close to a log-normal distribution characterized by one 
critical exponent ao = 2.28 ± 0.03 (the average over 130 critical states). Similar findings 
were presented already in Sec. 12.2 of Ref. 0. 



III. CORRELATIONS AT FIXED ENERGY 

In this section we first outline a theoretical description of scaling relations between the 
critical exponents of the correlations at a fixed energy and the f{ce) spectrum. We follow 
mainly the presentation of Ref. ||. Also, we present our numerical data which confirm the 
scaling relations. 

To study the spatial correlations of amplitudes for a fixed energy we consider the 
q— dependent correlations 

M^(r,L b ,L) := ([P l (L b )f[P l+s (L b )f) L (8) 

where the average is to be taken over all pairs of boxes with fixed distance r = sL b . 

For critical states where the microscopic scale (/# in our case) and the macroscopic scale 
(the localization length £ in our case) are separated, one can expect power law behavior of 
MM in the regime 

/ B <L b ,r,L<£. (9) 
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Usually, in critical phenomena one studies correlations for infinite system size (and Lb being 
microscopic) as a function of r alone. This is justified if, for large enough system sizes L, the 
correlation function is independent of L (for simplicity we neglect any trivial L-dependence 
due to prefactors in the definition of the observable P). However, this is not true in the 
multifractal case. Multifractality reflects broadness of the distribution function tt(P, L^/L) 
on all length scales. This is due to the dependence of the box-probability in each box on a 
large number of conditions, simultaneously. More generally, the local box observable Pi(Lb) 
depends on a large number of conditions for the entire system of linear size L, simultaneously. 
This behavior was denoted as "many parameter (MP) coherence" ||. In the context of the 
LD transition coherence at zero temperatures is due to quantum mechanical phase coherence 
of the electrons wave function, and disorder introduces a huge number of parameters, e.g. 
the position of point-scatterers. In the case of MP coherence one has to face the fact that 
MM depends non-trivially on L, even for L — > oo. To distinguish between the multifractal 
and the "ordinary" situations one can implicitly define a length scale L by the requirement 
that M^ will be independent of L for L > L. We call L the MP coherence length. Still, 
two different cases have to be distinguished. It may happen that L introduces a cutoff for 
correlations. For example, the correlation length £ is a MP coherence length of this kind. 
Alternatively, L does not introduce a cutoff and correlations still show a modified power law 
behavior with respect to r for r>L. Such kind of MP coherence length occurs in ordinary 
critical behavior where L is microscopic and no multifractality occurs, i.e. r{q) = d(q — 1). 

After these general considerations we can formulate our expectations for the present case 
of correlations of critical eigenstates at fixed energy. Because of the multifractal character 
of the critical states whose range is only limited by the localization length ( > i we can 
expect to be dealing with a MP coherent situation with L being the MP coherence length L, 
setting the cutoff for correlations. Therefore, we consider the regime lB^Lb<r<L<^t; 
and make the ansatz 

M'«l (r, L b , L) oc L x b 2{q) L- y ^r- z ^ . (10) 
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The task is now to find the scaling relations between the set of exponents £2(9), 2/2(9)5 z{q) 
and the r(q) function of Eq. (|3j). These scaling relations can be derived on the basis of 
heuristic arguments. Consider the limiting situations (i) where r is of the order of L b and 
(ii) where r is of the order of L. In case (i) the function will behave like ([P(L b )} 2( >) L 
while in case (ii) a decoupled behavior occurs ~ ({[P(L b )] q ) L ) 2 . The uniqueness of 
scaling exponents allows then to conclude the desired scaling relations 

y 2 (q) = d + r(2q) (11) 
x 2 (q) = 2d + 2r(q) (12) 
z(q) = d + 2r(q)-r(2q) (13) 

already proposed in p^J^j. It is worth mentioning that the sum £2(9) -2/2(9) — z{q) vanishes 
due to the normalization of the wave function. To verify numerically the scaling relations it 
is thus enough to verify two of the three equations given above. Therefore, we concentrate 
on the exponents z(q) and x 2 (q) for which we can set a fixed system size in numerical 
calculations. This reduces the computational effort substantially. 

Let us summarize the analytic behavior of x 2 (q) and z(q) according to Eqs. ( |T2"| , [IB]): x 2 (q) 
is a monotonic increasing function with negative curvature and asymptotic slopes given by 
2a(=Foo). It vanishes at q = 0. z(q) is non-negative with minimum at (0, 0). For q > 0(< 0) 
it is monotonically increasing (decreasing) and is asymptotically bounded by the dimension 
d (see Fig. |3|). To check on the validity of Eq. (|T3| ) we took 100 critical states of a system 
with L = 200l B and calculated M^(L b ,r, L) with fixed values L b = l B ,Al B ; L = 200l B . The 
distance r was varied from L b to 1501 B . The periodic boundary conditions in y-direction 
reduce the upper scale for a reliable determination of exponents to r <^ 100Z#- As can 
be seen from Fig. |] the power law behavior holds up to ~ 60/^- The numerical data for 
z(q) were obtained by determining the slope in the linear regime of the plots InM™ vs. 
lnr. In Fig. [5] the average of z(q) data over 100 states is shown in the regime \q\ < 2. For 
comparison the function d + 2r{q) — r(2q) is plotted, too. Within the errors the validity 
of the scaling relation Eq. (0) can be confirmed. For later comparison we mention that 
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z(l) = 0.43 ± 0.05, z(0.5) = 0.13 ± 0.03 and d 2 z/dq 2 (q = 0) = 1.1 ± 0.08 « 4(a - 0- To 
determine the exponents £2(9) numerically we fixed r = 100Z#, L = 200/# and varied L b 
from 21b to r. As is shown in Fig. ^|the scaling relation Eq. (|1^) is fulfilled, too. 

In summary, we have presented numerical evidence for the validity of the scaling relations 
Eqs. (PP). 



IV. CORRELATIONS AT TWO DIFFERENT ENERGIES 

In this section we firstly discuss the scaling relations for correlations of critical eigenstates 
at two different energies. We stress the role of a length scale related to the energy 
separation u. Secondly, we present numerical calculations which lead to the interpretation 
of L u being a cutoff scale for MP coherence. Thirdly, we verify the validity of scaling 
relations based on this interpretation. Finally, we discuss the relevance of our results for the 
interpretation of the local density of states as being an order parameter of the LD transition. 

We define the g-dependent correlation of box probabilities corresponding to two different 
eigenstates with energies E and E + uo 

MM(r,L b ,L) := ([Pi(E; L b )) q [P i+s (E + u; L b )] q ) L . (14) 

To understand the correlation behavior of non-localized states with respect to the energy 
separation one has to compare the relevant energy scales of the problem. These are the 
average level spacing A and the energy E c (u) corresponding to the time a wave packet 
(formed from states within an energy window of width to) needs to diffuse through the 
system, L 2 = (h/E c (u))D(u). Here D (uj) is the corresponding diffusion constant. According 
to Chalker 0, these scales give rise to the definition of two length scales depending on the 
energy separation uo: 

L w := (lo/E c (uj))- 1/2 L (15) 
^■^(u/Ar^L. (16) 



9 



The first length scale, L W) is the typical distance a wave packet will travel diffusively in 
a time U/uj. From this it is natural to assume that correlations between Pi(E;Lb) and 
Pi +S (E + u; Lb) will be present at least for distances r ^ whereas for larger distances 
the amplitudes are uncorrelated. Such uncorrelated behavior is typical for random matrix 
theory approaches to chaotic systems. The corresponding assumption is known as 'isotropy' 
or 'no-preferential basis' assumption and means that the unitary matrices that diagonalize 
the Hamiltonian are distributed uniformly in the unitary group and no correlations (apart 
from the unitarity property) between different matrix elements occur. Thus, the presence 
of correlations means a breakdown of the 'no-preferential basis' assumption. In electron 
systems with spatial disorder, however, a preference to some basis is always given. This 
preference is lost in M$ for distances r ^> L w . 

The second length scale, L u , is the linear size of a system with level spacing uo. Taking 
the existence of a preferential basis for being significant we will adopt the hypothesis that 
two wave functions with energetic separation smaller than the level spacing show a spatial 
correlation behavior of its amplitudes similar to that corresponding to one of those wave 
functions, i.e. they are statistically indistinguishable. 

At the critical point of the LD transition the conductance becomes a size independent 
quantity, G = g*e 2 /h, and with the help of the Einstein relation between conductivity and 
diffusion one finds L w = (g*) l l d L w |§. Since g* is of 0(1) the two length scales coincide at 
the LD transition. Therefore, we will focus our considerations on the role of L u and start 
with the following working hypothesis: 

M$ will show the same correlation exponent z(q) as Mj%L provided r < L u . For dis- 
tances r ^> the amplitudes Pi(E; L b ) and P i+S (E + u>; L b ) are uncorrelated. 

This means that is a MP phase coherence length setting a cutoff for correlations. Asking 
for the scaling properties of M$ in the regime L b < r < L w < L we make the ansatz 
(cf. Eq. (0)) 

MM oc L ^) r -^) L H q ) L -Y 2 ( g ) (17) 
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Here we have already anticipated that the exponent with respect to r is z(q), as given before 
and will check if this is consistent within the following procedure of deriving the scaling 
relations between X 2 (q), z(q), Z(q), Y 2 (q) and r(q). As in the case of for a fixed energy 
we consider limiting situations. In case (i) we keep r of the order of Lb and of the order 
of L, whereas in case (ii) we keep r of the order of L^. Following our working hypothesis 
for case (i) we expect no difference to the case with zero energy separation 

M^<x([P i (E;L b )]^) T , (18) 



/ L 

resulting in 

X 2 (q) - z(q) = d + r(2q) = Y 2 (q) - Z(q) . (19) 

For case (ii), according to the hypothesis, M$ depends neither on r nor on and ap- 
proaches the uncorrected value, 

M^^m{E;L b )f) L f . (20) 

This leads to 

z(q) = Z(q) ; X 2 {q) = 2d + 2r(q) = Y 2 {q) . (21) 

Equations ( |19| , pTD yield the scaling relations 

X 2 (q) = 2d + 2r(q)=Y 2 (q) (22) 
z(q) = Z(q) = d + 2r(q)-r(2q) (23) 

which form a central result of this article. Now the following conclusions can be drawn: 
1. The result for z(q) is the same as in the case of zero energy separation. 2. The energy 
separation u is not an independent scaling parameter but appears only in the combination 
L^/r. 3. The exponent corresponding to the box size, x 2 (q) = X 2 (q), is not affected by 
a finite energy separation but the exponent corresponding to the system size L (which is 
y 2 (q) for zero energy separation) splits up into the exponents Z(q) (corresponding to LJ) 
and Y 2 (q) (corresponding to L for finite energy separation). 
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In Ref. it has been speculated that for r ^> L w modified power law correlations may 
still exist. Such behavior would be in conflict with the interpretation of adopted here 
and it would lead to different scaling relations in comparison to Eqs. ( |2~2"| , |2"3| ) . Especially, 



the equality between z(q) and Z(q) would be lost for arbitrary values of q. However, the 
possibility of such behavior could not be ruled out on the basis of previously obtained 
numerical data which correspond to the case q = 1. 

We will now demonstrate that the working hypothesis formulated above and the resulting 
conclusions 1. and 2. are consistent with our data. 

Firstly, we checked if serves as a cutoff for correlations. That this is indeed the case 
can be seen from Fig. [7] showing the logarithm of M$ as a function of In r (r is measured 
in units of Is) for a fixed size L& = Ib and q = 0.5. The energy separation u corresponds 
to = 301b- A clear r-dependence up to approximately this scale can be observed. In 
order to calculate a reliable scaling exponent z(q) one has to take > 1001b- Doing so, we 
found (within the errors) the same z(q) values as for correlations with no energy separation, 
e.g. z(l) = 0.39 ± 0.04, z(q = 0.5) = 0.10 ± 0.03. Thus, the working hypothesis turned out 
to be consistent with our data. The next step is to investigate the scaling exponent Z(q) 
corresponding to which should, according to the scaling relation Eq. fl2"3]), be equal to 
z(q), resulting in a combined scaling parameter r/L w . 

For the numerical determination of Z(q) we used about 130 2 combinations of pairings 
between different eigenstates. For each pair the actual value of the cutoff scale for correla- 
tions fluctuates around the calculated value L^. This fact requires a large amount of data to 
extract reliable scaling exponents. From the plots of In M$ versus In L w with fixed values of 
r « Lb = Ib we determined the approximate linear behavior in a regime between « 201b 
and L w « 501b- As can be seen in Fig. [?| (q = 0.5) there are fluctuations due to rare pairings. 
In Fig. [8] we show the numerically obtained Z(q) function together with the errors of the 
linear fit in log-log plots. It is compared with the previously obtained z(q) function for zero 
energy separation (cf. Sec. [mp . The data show that the equality Z(q) = z(q) is consistent 
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with our data, e.g. Z(l) = 0.38 ±0.04, Z(0.5) = 0.09 ±0.03 and d 2 Z/dq 2 (q = 0) = 1.0 ±0.1. 

We have thus demonstrated that the working hypothesis as well as the scaling relations 
are consistent with our numerical results. 

Having established the role of r / L w as the relevant scaling parameter for correlations of 
eigenstate amplitudes with universal exponent z(q) (related to r(g) by a scaling relation) 
let us now discuss the consequences of these findings for the interpretation of the local 
density of states being an order parameter of the LD transition. The local density of states 
(LDOS) is defined formally as p(E, r) = J2n 3(E—E n ) \ip n (r)\ 2 where ij} n {v) are wave functions 
corresponding to eigenenergies E n . In a finite system this function has isolated peaks at E n 
and becomes a smooth function of energy only for an open system or in the thermodynamic 
limit. This is also true for the global density of states (DOS) p(E) = L~ d J2 n 3(.E ~ En)- 
Since the smallest relevant energy scale for the structure of the DOS as well as of the LDOS 
is set by the average level spacing A a smearing out of the (^-functions over this scale is 
needed to talk about DOS and LDOS in finite systems. It is known that the average DOS 
does not reflect the LD transition but is a smooth function of energy and independent of 
system size L. With this smearing out of the 5-functions we define the LDOS as 

p(£,r) = A(£)-V(£,r)| 2 (24) 

where \i/j(E, r)| 2 stands for the microcanonical average of squared amplitudes at a given 
energy E. Since A(E) is smooth in energy and behaves as L~ d , the scaling behavior of the 
LDOS is determined by that of the wave function. Consequently, we have 

([p(E,v)Y) L ^L^ d -^ (25) 

and for the typical value 

p typ = exp[(ln(p(r))) L ] oc L d ^ (26) 

which does reflect the LD transition. Scaling L with the localization length £ oc \E — E*\~ u 
where E* and v are the critical energy and the critical exponent of the localization length, 
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respectively, we arrive at the conclusion that the typical LDOS vanishes on approaching the 
critical energy with exponent /3 typ = u(a — d). This observation has led to the interpretation 
of the LDOS being an order parameter field of the LD transition [||] (see also P,|T5f). The 
unconventional feature as compared to ordinary critical phenomena lies in the facts that (i) 
the order parameter field has a broad distribution resulting in a non-linear dependence of 
exponents on the degree of moments considered (multifractality) and (ii) that the average 
value has vanishing scaling exponent while the typical value gives rise to a positive scaling 
exponent /3t yp . The scaling relations that we derived for the wave functions amplitudes 
transform to scaling relations of the LDOS since each box amplitude has to be multiplied 
by a constant factor of L d , 

((p(E, v^YipiE + u, v 2 )Y) L cc (Ljr)^L- A ^ ; r = |n - r 2 | (27) 
z{q) = d + 2r{q) - r{2q) ; A{q) = 2(1 - q)d + 2r(q) . (28) 

In this light our scaling relations turn out to be the appropriate scaling relations connect- 
ing the spatial correlations of the local order parameter field to its scaling dimensions (cf. 
Eqs. (|2B]),(P7|),([28])). However, due to the MP coherence these scaling relations are different 
from those of ordinary critical exponents where only one correlation exponent appears and 
no multifractality has to be taken into account. The most remarkable difference lies here in 
the observation that power law scaling is also present in the MP coherence length L. The 
latter is given by the system size L for energy separation smaller than the level spacing 
or by L w for energy separation larger than the level spacing. Our findings demonstrate 
that the interpretation of the LDOS being an order parameter field for the LD transition is 
supported by the existence of universal scaling relations. We close this section by pointing 
out that z(l) = 2 — r(2) « 0.4 ^ (for the correlator of the density of states) with L^/r 
forming the scaling parameter is equivalent (cf. PJIBfl) to the phenomenon of "anomalous 
diffusion" found by Chalker and Daniell || . As pointed out in |8j , the anomalous character 
of diffusion lies in the non-Gaussian dispersion of a wave packet in time t despite the fact 
that the average diameter grows like y/i. This non-Gaussian time dispersion is caused by 
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the multifractal character of eigenstates. 



V. CONCLUSIONS 

We have demonstrated explicitly that the distribution of amplitudes of critical eigenstates 
in quantum Hall systems (QHS) is contained in the multifractal f(a) spectrum, essentially 
characterized by one critical exponent, «o ~ 2.28 (Sec. |T|). Following Ref. || we derived 
scaling relations which relate the critical exponents of the spatial correlation of amplitudes 
(for a fixed energy in the critical regime) to the multifractal spectrum (Eqs. (p7T| , [T~2| , p~3D ) and 
demonstrated that they are consistent with numerical data obtained for QHS (Sec. 
Most interesting is the (non-trivial) dependence of the correlator on the system size. In 



section [TV] we considered correlations of amplitudes corresponding to critical eigenstates 
with energy separation u. Following Chalker || we identified a length scale (describing 
the system size with level spacing u) as the relevant cutoff for correlations. Furthermore, 
we exploited the hypothesis that the amplitudes are correlated as for zero energy separation 
provided their distance is much less than L w . We found scaling relations (Eqs. (p5|, f2T| , f28| )) 
which relate all the correlation exponents to the universal multifractal spectrum of critical 
eigenstates. Most important are (i) the confirmation of the scaling parameter L^/r and 
(ii) the identification of as the upper limit for multifractal correlations. We discussed 
implications of our results for the statistical properties of the local density of states in the 
critical regime. Our findings demonstrate that the interpretation of the LDOS being an 
order parameter field for the LD transition is supported by the existence of universal scaling 
relations. 
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FIGURES 

FIG. 1. Squared amplitudes of a critical wave function for a system of linear size 200 magnetic 
lengths. The corresponding f{a) spectrum (•) together with the parabolic approximation (• • •) 
(PA, Eq. (^)) and the semi-elliptic approximation (-) (SEA, Eq. ([?])) are also shown. 

FIG. 2. Histogram of the logarithm of squared amplitudes shown in Fig. [I]. The continuous 
curve is the distribution function following from the corresponding f(a) spectrum via Eq. (||). 

FIG. 3. The function of critical exponents z(q) following from the the scaling relation Eq. (|l5|). 

FIG. 4. The correlator AfM (for q = 0.5) as a function of the spatial distance r in a log-log plot. 

FIG. 5. The numerically obtained data for z(q) (*) in comparison with data following from the 
scaling relation Eq. ( |i~3| ) (o). 

FIG. 6. The numerically obtained data for X2(q)- The line shows the data following from the 
scaling relation Eq. (|i~2|). 

FIG. 7. The correlator AlL q] (for q = 0.5) function of the spatial distance r in a log-log 
plot for = 301b- 

FIG. 8. The correlator AfJ? 1 (for q = 0.5) as a function of the length scale L u in a log- log plot. 

FIG. 9. The numerically obtained data for Z(q) (o) in comparison with the numerically ob- 
tained data for z(q) (*). 
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